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We examine four WTC 7 descent curves, labeled "C," "E," "N," and "O," ei- 
ther anonymously published, or confidentially communicated to us. Descent curve 

o ■ 

CN ' describes apparent height of a collapsing building as a function of time. While all 

sets are mutually consistent, it is set "C" which suggests that there are three active 

hj ; 

phases of collapse. Phase I is a free fall for the first Hi ~ 28 m or T\ — 2.3 s, during 
which the acceleration a is that of the gravity, a = g = 9.8 m/s 2 . In Phase II, which 
continues until drop Hi ~ 68 m, or T2 — 3.8 s, the acceleration is a ~ 5 m/s 2 , while 

1 ' 

O in Phase III which continues for the remaining of the data set, a ~ — 1 m/s 2 . 

o : 

Qh, We propose that the collapse of WTC 7 is initiated by a total and sudden anni- 

t/5 ■ 

hilation of the base (section of the building from the ground level to -Hi), which 

• i-H ' 

C*2 ■ 

£>y then allows the top section (building above Hi) to free fall during Phase I, and 

Q_i, then collide with the ground in Phase II and III. We interpret the latter two phases 

^j. I of the collapse as the top section being comprised of two zones, the 60% damaged 
> 

■ primary zone (below H2) and the intact secondary zone (above H2). We derive a 

G\ 

t^. . physical model for collision of the building with the ground, in which we correct 

i the "crush-up" model of Bazant and Verdure, J. Engr. Mech. ASCE, 133 (2006) 

qq ' 308. The magnitude of resistive force in the two zones of the top section obviates 

O ■ 

the catastrophic failure mechanisms of Bazant and Verdure (ibid.), and of Seffen, J. 
Engr. Mech. ASCE, 134 (2008) 125. The total duration of the collapse, assuming 
I that Phase III continues to the end, is in the range 7.8 — 8.6 s. 

We compare our findings to those of NIST investigators and find an agreement with 
respect to the distribution of damage in the primary zone. We conclude that the 
building was destroyed in a highly controlled fashion. 
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I. INTRODUCTION 
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FIG. 1: Position of the top of WTC 7 as a function of time (descent curve) as recorded in four 
data sets, labeled "C," "E," "N," and "O." Data set "O" covers some 1.8 s of collapse, however, 
it extends to 9 seconds before the collapse and so shows that during that time the building did 
not noticeably move. Data sets "E" and "N" cover the first 3.5 seconds of collapse with different 
sampling rates. Most of the report deals with "C" as it is the most extensive: it covers 4.8 seconds 
of collapse. 



World Trade Center (WTC) 7 perished after WTC 1 and 2 collapsed on September 11, 
2001. Its demise was examined by Federal Emergency Management Agency (FEMA), which 
in 2002 issued a report, jlj In the report it was claimed that the most likely cause of the 
collapse was the gradual weakening of vertical columns in the lower part of the building 
following their exposure to the mechanical and thermal stress. As the sources of mechanical 
stress the falling debris and the earthquake from the collapse of WTC 1 and 2 were given, 
while the fires raging inside the building, some of which were fueled by the heating oil known 
to have been stored in the building, were cited as a source of thermal stress. In November 
20, 2008 the NIST investigators issued a report in which it was claimed that the fires that 
followed the impact of debris from the collapse of WTC 1 (the north tower) led to the 
collapse of WTC 7. Q 

One easy-to-capture feature of a collapse of any high rise building is the motion of its 
top as a function of time - its descent curve. If the falling building collides with the Earth 
at a relatively well defined collision plane, then, on one hand side, the motion of the top 
reflects the motion of the moving part of the building. On the other, this feature allows us 
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to write down the equation of motion for the moving part. The physical model that follows 
from the equation of motion, can be fitted to the descent curve, which in turn allows us to 
estimate how the collapse was initiated and when, and what was the damage distribution in 
the building. 

The goal of this report is the analysis of motion of the top of WTC 7 as recorded in four 
descent curves: "C"|a|, "E"fl, "N" 5|, and "O" 6], cf. Fig.HJ We develop a physical model 
of the descent and fit it to the curves in order to obtain a local building's strength. We 
compare so obtained values to the estimates provided by us [7j, by Bazant and Verdure [8j 
and by Seffen We propose a scenario of collapse initiation which is consistent with the 
descent curve in its entirety, and which distribution of damage prior to collapse agrees with 
that of the NIST investigators. ^ 



II. DESCENT CURVE AND ITS PHYSICAL MODEL 
A. Finite Differences Analysis of Descent Curve "N" 



TABLE I: Descent curve "N"|a| and its acceleration as found by using the finite differences. 
Shown are the results for the relevant velocities: the mean velocities on the intervals (vi for interval 
and the momentary velocities at the ends of the intervals, Vi = v{tj). As can be seen the 
average acceleration Oj for the first T± ~ 2.3 s, or the displacement of Hi ~ 26 m, oscillates around 
the gravity, g = 9.81 m/s 2 , indicating a free fall. After time T\ the acceleration drops to ~ 5 m/s 2 . 
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The descent curve shows the motion of the top of the building as a function of time. We 
use the term "top section" to label the part of the building which motion is common to 
that described by the descent curve. Our analysis of descent of WTC 7 then has two goals: 
to find the extent of the top section, and to determine its acceleration. We recall that the 
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acceleration of the top section is a result of the forces acting on it. We discuss the forces 
acting on the top section later when we derive a physical model of the event. 

We start with data set "N" because it is short: it contains 8 points at 0.5 s apart. A 
reader interested in full presentation of how the data was obtained is kindly directed to the 
original publication [5J]. We note that the estimated error of the distance is ±0.5 m. 

The results of the finite differences analysis of the descent curve "N" are given in Tbl.[T] 
together with two intermediate velocities: the mean velocities Vi on intervals and 
the momentary velocities at the end points tj. We immediately notice that the mean 
acceleration of the top section has two distinctive values, which we associate with the phases 
of descent and label with Roman numerals. During Phase I the average acceleration is 
a = 9.69 m/s 2 , which is within 1% from the free fall acceleration given by the gravity 
g = 9.81 m/s 2 . This phase lasts for the first 2.35 — 2.85 s. Phase I is thus a free fall phase. 
Phase II continues during which the mean acceleration, a, drops to a ~ 5.1 m/s 2 . This 
phase presumably continues for the remaining one second of recorded data. 

Now that we have determined the acceleration of the top section, we answer the question 
of how far below the top of the building does the top section extend. We notice that the 
2.35 s long free fall corresponds to a distance of ~ 26 m. We label this 26 m section of 
the building in the path of the top section "the base." The question thus becomes at what 
height is the bottom of "the base." The clue about where is "the base" comes from the 
NIST investigation on the collapse of WTC 1 and 2, where it is explicitely hypothesized [h]] 
(p.146, Sec. 6.14.4) 

The structure below the level of collapse initiation offered minimal resistance to 
the falling building mass at and above the impact zone. The potential energy 
released by the downward movement of the large building mass far exceeded the 
capacity of the intact structure below to absorb that through energy of deforma- 
tion. 

In other words, when a steel frame (the top section) collides with another steel frame ( "the 
base") it feels minimal to no resistance. However, as indicated by the descent curve, when 
the top section finishes destroying "the base" its acceleration suddenly drops to half its 
value. We conclude that at that point the top section collides with the only other object a 
falling body can collide with - the Earth's surface. This puts "the base" to the base of the 
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building, from the ground level to the height Hi ~ 26 m, and the top section from Hi all 
the way to the top. This positioning of "the base" is consistent with the visual appearance 
of the collapse, where the visible part of the building moves uniformly. Behavior of the base 
of the building, on the other hand, is completely hidden in the cloud of dust created at the 
onset of collapse. 

Here however, the following comment is due. The visual appearance of collapse, and 
positioning of "the base" to the base of the building do not imply that the NIST hypothesis 
stated above is correct. Had the top section started its motion from Hi ~ 26 m down and 
accreted the mass of the base in its path this would create an inertial brake which would 
prevent the top section from ever achieving the free fall acceleration. As the top section 
starts immediately with the free fall acceleration this implies that (i), the base is converted 
into a free falling rubble prior to the top section falling through it, and (ii), the crushing of 
the base is not done by the top section. We return to this point later in the report, when 
we discuss it in more quantitative fashion. 

Based on this argument we assert that the collapse of WTC 7 begins with a sudden and 
total annihilation of the base, which, on one hand, allows the top section to free fall for the 
height of the base, and on the other, leads to the observed change in acceleration once the 
top section reaches the ground. We label euphemistically the collapse initiation moment a 
"release," and the height Hi ~ 26 m from which it occurs a "release point." While Phase I 
describes the free fall of the top section to the ground following its release, in Phase II the 
top section collides with the ground. Bazant and Verdure [8| label Phase II of the collapse 
a "crush-up," and we adhere to their terminology for the rest of the report. 



B. Physical Model of Collision of a Building with the Ground 

1. Derivation of General Model 

Use of an one-dimensional model for description of a collapsing building is justified pro- 
viding that the building of interest, WTC 7, collapsed almost perfectly in its footprint. This 
remarkable feature allows us to exclude any transverse coordinates from the analysis and 
use only the height to describe the motion of the building. In doing so, we in effect average 
the behavior of the building with respect to the excluded coordinates. Most importantly, we 
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FIG. 2: "Crush-up" model of collapse. The top section, between drops Z and Z%, of the building 
collides with the ground at a well defined collision plane at Z\. In the "crush-up" the mass and 
the momentum of the top section are lost by being transferred to the ground. Finite linear density 
of the crushed building, p > 0, leads to a slow upward crawl of the collision plane so that Z\ is not 
exactly at the ground level. A resistance the top section offers to its destruction at the collision 
plane is described by a resistive, or "crushing," force. 

do not need to consider a failure of individual load-bearing structural elements, e.g., vertical 
columns, and instead concentrate on their collective response. 

The collapse dynamics of the building in "crush-up" mode is shown in Fig.[2j The falling 
building, the top of which is at Z, collides with the ground at a collision plane at Z\. As 
a result of a finite compaction ratio the collision plane crawls up, toward the top section. 
First, we assume that the building is of uniform mass density po = M/H, where M is the 
total mass of the building and H its height. We introduce a compaction ratio k, as 



where p is the density of compacted building. For simplicity, we assume that the compaction 
is uniform as well, i.e., that p is a constant. Second, we introduce two coordinates to mark 
the progression of "crush=up," an apparent drop of the top of the building, Z, and, a 
position of the collision plane, Z\. The two coordinates are connected by the requirement 
that the mass of the building is conserved, 




(2) 
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yielding 



and 



Z 1 = H- j^—Z, (3a) 
Z 1 -Z = H- -^—Z, (3b) 

1 — K 



z > - "TV" (4) 



We proceed with the derivation of equation of motion for the apparent drop of the building 
Z = Z(t). This is most easily accomplished by using the energy formalism. The kinetic 
energy of the moving part of the building is 

K(Z, Z)= l -p, (Zx - Z) Z 2 = \ po (H - - J— Z) Z\ (5) 
z z 1 — /c 

while its momentum is 

The potential energy of the building is given by U (Z) = — f!? dX p(X) g X, giving 

U{Z) = - l -pg {H 2 -Zl)- l -p Q g (Z 2 - Z 2 ) = -^p g (h 2 + 2 H Z - -^j. (7) 

The gravitational force with respect to the coordinate Z follows from G = —dU / dZ, and is 
given by 

G = p -g-(H- T ^-Zy (8) 

Last energy in the problem is the "latent" energy L, from which the force with which the 
building resists its destruction, call it resistive force R, is derived. We have, 

L = -[ 1 dXR{X) = L{Z X - Z). (9) 
Jo 

The resistive force R with respect to the coordinate Z follows, as before, R = —dL/dZ, and 
is given by 

= _az,(z, - z) = - z) S(Zl - z) = _j_ / _ _._ \ 
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This said, the equation of motion for Z follows from Newton's law, 

P = G + R+(P) . (11) 

V / loss 

The loss of momentum (mass, energy) occurs at the avalanche front where the momentum 
is transferred to the stationary part of the building. The loss rate is Z ■ m, where m is the 
mass of the moving part, yielding for the equation of motion, 

1 R(H-^-Z) 
Z = g + — ) 1 -=?-L. (12) 

We observe that while in the limit k — > Eq. (1121) coincides with the result of Bazant and 
Verdure 8j, for k / their model does not correctly incorporates compaction. 

The resistive force R describes how the building resists its destruction at the avalanche 
front. It is a function of strength of the structural elements of the building, as well as their 
failure mode. Most notable contribution comes from the vertical columns, the strength of 
which varies with height Z. For simplicity, we assume that the dependence of R on Z is at 
most linear, yielding 

«<*>-,. (r + 4l (13) 



where r and s are two dimensionless parameters. With this parameterization of R we obtain 
the ordinary differential equation (ODE) for Z, 



* = °-V-T^- a-Ja T * *» • (14) 



1-K H' 



Finally, we note that if we set r = s = 0, we immediately obtain the equation of motion for 
a free fall, 

Z = g. (15) 



It has to be kept in mind that Eq. ( 1141) contains an assumption of what happened to the 
base at the onset of collapse: it was instantaneously converted into a free floating pile of 
debris, which for t > starts to free fall to the ground and contributes to motion of the 
collision plane at Z\. 

When we solve Eq. f fMl) later, we always assume that the top section starts its motion 
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from rest, Z(0) = Z(0) = 0. 

Parameters of the building that enter Eq. ( fT4l) are the total height of the building, H = 
186 m, and the parameters r and s of the local resistive force, R/(M g). Furthermore, in 
our simplified model the building comprises 47 floors, each 3.66 m high, and the 14 m high 
lobby. 

Though it might not be obvious, the compaction parameter k is of secondary importance. 
Thus, in what follows we further simplify Eq. (j!4p by taking k = 0. 



2. Apparent Weight of the Building During the Collapse 

An apparent weight the top section exerts during the collapse on the Earth's crust, 
W'/(M g), is given by 



W 
~Mg~ 



%, for t < 0, 



forte [0,71), (16) 
f + |i 2 + fa + Sj ■ (1 - z)) for t>T x . 



where z = Z/H is the scaled drop, and z = Z/(H/T ) is the scaled velocity, with T 2 = 2 Hj g 
being a free fall time from height H . Index j in rj and Sj keeps track of the phase of collapse 
so it is a function of the time, as well. 

In Eq. (1161) . M'/M = 1 — H\j H is the mass of the top section, T\ ~ 2.3 s is the duration 
Phase I, the free fall, while Hi ~ 26 m is the distance. The terms appearing in W'/(M g) 
are, from left to right, (z), the weight of the top section that has already reached the 
ground, pg(Z — Hi); (ii), reaction force due to the change of momentum of the crushed 
material at the collision plane, pZ 2 ; and (Hi), the resistive force at the collision plane, 
R = R(H — Z). The last term is present because the crushing of the building at the collision 
plane is performed between the ground and the top section. 

We believe that W, and in particular its time derivative, can be used in interpretation of 
the seismic signal of the building's collapse. As an attempt to connect the two brings forth 
numerous additional complications which need to be properly addressed, we leave this topic 
to future publications. 
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III. DESCENT CURVE "C 




0123456789 10 



Time / seconds 

FIG. 3: Results of fitting to data set "C" (black points) to a three-zone model (position in red, 
acceleration in orange). Descent starts with a free fall for a distance of Hi ~ 28 m and lasts 
T\ = 2.35 s. In ensuing collision with the ground varying acceleration reveals that the top section 
consists of two zones, the primary below Hi — 68 m, for which a ~ 6 m/s 2 , and the secondary 
above, for which a ~ — 1 m/s 2 . 

Data set "C," as presented to us, consists of two descent curves, one at 10 samples per 
second, and the other at 5 samples per second. We leave presentation and discussion of the 
descent curve to future publication by their author. [3] In what follows we use low resolution 
data set only, which consists of 25 data points, with an estimated error in distance of ±0.2 m. 

The importance of the physical model ( fl4l) is that it allows us to identify the phases of 
collapse (stages of descent between which the acceleration changes discontinuously) to the 
zones of the top section being destroyed in the collision with the ground. 

As a measure of how close is a solution of Eq. (I14p . call it Z = Z(t), to the data set "C" 
we use a sum-of-absolute-errors (SAE), 

N 

SAE=J2\Z(U;{r J ,s J } J )-z l \. (17) 

i=l 

Inspection of data suggests that there are three zones in the building, each with its own r 
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and s. 

We estimate {(r i; Si)} i= i j3 as follows. We perform a sequential fitting where in the first 
step we fix heights Hi and H 2 while varying simultaneously the parameters of the resistive 
force in all three zones. In the second step we vary the positions of one of the bounds, Hi 
or H2, in 0.5 and 1 m increments, respectively, where we use the results obtained in the 
previous step as the initial conditions. For minimization of SAE we use a simplex method 
of Nelder and Mead, as it does not require computation of derivatives. For solving the ODE 
we use a Runge-Kutta Prince-Dormand (8,9) method. Both methods are implemented in 



the GSL. 
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The results of the minimization are shown in Fig.[3j In the three-zone model the best 
fit achieves SAE ~ 3.824 m. The boundaries between the zones are at Hi = 28 m and 
H 2 = 68 m, where (r 1; si) = (0,0), (r 2 ,s 2 ) = (0,0.41) and (r 3 , s 3 ) = (0,1.09), for Phases 
I, II, and III of descent, respectively. The fit of the model to the data is excellent with an 
average error per point being ~ 0.15 m, thus smaller then the ~ 0.2 m margin given by 
the set "C" author. We estimate the width of the boundaries by introducing a number of 
"micro" zones at Hi and H 2 . This procedure yields for the width of transition at Hi of 1 m 
and 2 m at H 2 , which we write as Hi = 28 ± 0.5 m and H 2 = 68 ± 1 m. The accelerations 
during the descent are a ~ g for the free fall phase (Phase I), a ~ 6 m/s 2 for the "crush- up" 
of the primary zone (Phase II) and a ~ —1 m/s 2 for the "crush-up" of the secondary zone 
(Phase III). 

We estimate uncertainties in {(rj, Si)}i=i^ as follows. We perform a large number of 
optimizations (N = 100), where we randomly choose initial values for (rj, Sj), while keeping 
Hi and H 2 fixed. Further, of all so obtained {(rj, Si)}i=i^ we keep only those which SAE is 
within 1% of the best (smallest) value. In Fig. H] we show the results of minimization. While 
(r 1; si) remain unremarkably close to 0, in Phase II we find that (r 2 ,s 2 ) , s spread along the 
line, 

r2 +-^1. (18) 



0.31 0.41 

For Phase III we find that (r^, ss) y s spread along the line 



+7^1. (19) 



0.61 1.10 



This indicates that SAE posses shallow minima along the lines Eq. ( |T8|) and Eq. (|T9l . 
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FIG. 4: Uncertainties in (r, s) in the three-zone model obtained by allowing SAE, Eq. (I17j) . to vary 
up to 1% from the minimum. As a result (r^, Sj), i = 2, 3, disperse along two straight lines, Eq. <\18\) 
for Phase II, and Eq. (I19p for Phase III. At the same time, the duration of collapse spreads uniformly 
between 7.5 and 8.4 seconds. For comparison we also show our estimates for r and s in WTC 1 
and 2 (black circles) from jjj: intact building with (r, s) = (0.2,0.7), and the contribution from its 
perimeter columns which (r, s) = (0.05,0.3). The radii represent their uncertainties: ~ 20% for r, 
and half that for s. We see that the secondary zone of WTC 7 appears to be intact. 

We further reduce a number of zone parameters in the three- zone model by considering 
that the top section throughout both zones is comprised of load bearing elements of identical 
properties. Then, the anticipated difference in strength between the zones comes from 
varying their number of elements. We introduce a constant k, which couples the resistive 
force in the primary and the secondary zone, or in Phase II and III, respectively, as follows 



While the results of optimization remain unremarkably close to the one listed previously, 
this procedure yields a very narrow estimate k = 0.42 ± 0.04. In other words, the primary 
zone appears to be 60%, or so, compromised compared to the secondary zone. 




(20) 
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A. The Secondary Zone 



To our knowledge, there were no attempts to estimate the magnitude of resistive 



WTC 7. However, such estimates were provided for WTC 2 by Bazant and Verdure 



orce m 
, and 

for WTC 1 and 2 by us [7]. The following discussion is based on an assumption that those 
estimates represent reasonable values for WTC 7, as well: 



Bazant and Verdure's estimate: The authors make an educated guess that the 



crushing energy AL per floor in an intact building is AL = 0.6 — 2.4 GNm. [14[ The 
resistive force is then given by R = AL/AH, where AH = 3.7 m is the floor height. 
They take mass of WTC 1 and 2 to be M = 3.2 • 10 8 kg which yields r = R/(M g) = 
0.05—0.2. The authors defend their estimate being so small by introducing a scenario of 
collapse which is currently being disputed. Seffen Q proposes an alternative, more 
obscure, catastrophic mechanism, which net effect is the same: a near-zero resistive 
force in a collapse of otherwise an intact building. 

Given the magnitude of the resistive force in three phases of collapse, the two catas- 
trophic mechanisms can, at best, be applied to the rapid reduction of the base's 
strength to zero in Phase I. We note that the physical processes behind both mech- 
anisms require the top section to produce "destruction waves" at the collision plane 
with the base. The "waves" propagate through the structure in front of the top section 
and reduce its resistive force to near-zero, after which the top section collects the pieces 
in its path. However, this cannot explain a free fall seen in WTC 7: the collected mass 
acts as an inertial brake which prevents the top section from ever reaching the free fall 
acceleration. We quantify this analysis in the next section. 

We conclude that their physical mechanisms, and consequently their estimates of re- 
sistive force, do not apply to WTC 7. Whether they apply to WTC 1 and 2 is yet to 
be seen. 

Beck's estimate: we argued that it suffices to consider a textbook model of resistive 
force, 16j R = e ■ Y , where Y = Y(Z) is the ultimate yield strength of the vertical 
columns as a function of drop from the building's top, while e ~ 0.25 is the ultimate 
yield strain of the structural steel used in the building. Using the estimates for prop- 
erties of the vertical columns in WTC 1 and 2 (their cross section, strength and the 
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weight they carry) one arrives to Y/(M-g) ~ 0.8 + 2.7- Z/H, that is, (r, s) ~ (0.2,0.7) 
in an intact building. On the other hand, the parameters of the secondary zone of 
WTC 7 are on the line + S3/I.I — 1. We see that, within a margin of error 

of under 5%, the two overlap. We base our error estimate on the ratio of a distance 
between the two (~ 0.04) to r + s = 0.9, and which is under 0.05. This is illustrated 
in Fig.H] where we plot r and s in the primary and the secondary zone in WTC 7 
and compare it to their values in intact WTC 1 and 2, and to the contribution from 
perimeter columns only. 

We conclude that the secondary zone of WTC 7 appears to be intact. 

B. The Primary Zone and the NIST Hypothesis Regarding the Collapse Initiation 

The descent curve provides us with a precise estimate of the primary zone being 60%, 
or so, damaged compared to the intact secondary zone. Here, the primary zone stretches 
from Hi = 28 m (4th floor) to H = 68 m (15th floor). We use the findings by the NIST 
investigators to posit that the core columns are absent in the primary zone. 

We start by stating the hypothesis of the NIST investigators, [2] regarding the distribution 
of damage prior to the collapse and its initiation, both illustrated in Fig.0 The report 
blames a failure of the column No. 79 near the 13th floor for the collapse initiation. Its 
failure presumably induces a cascade of floor failures, which cause the buckling of "additional 
columns" and "within seconds, the entire building core is failing." 

We immediately note that the last statement in their hypothetical failure scenario is not 
correct: according to the descent curve it is not the entire building core that is failing but 
only the core below H 2 . We recall that the descent curve indicates that the secondary zone 
(the building above H 2 ) managed to stay intact not just at the initiation of collapse but also 
until the collision plane of destruction reached it some 4 seconds into the collapse. Now, if it 
is only the core below H 2 that is failing, than this failure for sure includes the core columns 
of the primary zone, and possibly the core columns of the base. This is only the first of 
disagreements between the NIST hypothesis and the descent curve. 

Second, it is not clear what drives the failure of the core below H 2 . It is not the weight of 
the building above: the core columns are severed at H 2 so below H 2 they carry little more 
then their own weight. Thus, the mass participating in the "cascade of floor failures" at 
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FIG. 5: Status of the WTC 7 primary zone at the onset of collapse as proposed by the NIST 
investigators. The failure of column No. 79 near the 13th floor is presumably responsible for the 
initiation of the collapse. 0] 

best corresponds to the mass of the core below H 2 . 

Third, a connection between the cascade of floor failures and the release of the top section 
at H\ ~ 28 m is not clear in the findings by the NIST investigators. The hypothetical failure 
that starts near H 2 propagates downwards through the center, possibly to the ground level, 
on one hand side. On the other, at H 1 this "cascade" surfaces at the sides of the building 
and severs the perimeter columns so that the free fall may commence. If the "cascade" 
indeed surfaces near Hi, and is "spontaneous," then Hi should vary considerably at the exit 
points along the perimeter of the building. In terms of the descent curve this would manifest 
itself as Hi being a range over which the resistive force changes from zero to some other 
value. Contrary to being "spontaneous," the descent curve indicates a sudden transition 
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from Phase I to Phase II, where its width is less than 1 m. 

We conclude that the NIST report, as is, agrees with the descent curve only in regard 
to the distribution of damage in the primary zone: in a cascading floor failure that started 
at the top of the primary zone the building, most likely, lost all of its core columns in the 
primary zone and in the base. We emphasize this point in Fig. H] which, among others, shows 
that r and s in the primary zone of WTC 7 are rather close to the estimated contribution 



from the perimeter columns in WTC 1 and 2, we presented in jj]. 

The NIST investigators, on the other hand, ignore the demise of the base and a role it 
played in initiation of the collapse. 

C. Free Fall and the Demise of the Base 

From the beginning it transpired that the top section flew through the base in a fashion 
awfully close to a free fall. Later when discussing the secondary zone we mentioned an 
alternative physical model of descent: an actual avalanche. We now examine more closely 
Phase I of the descent for a distance Hi ~ 28 m as described by the data set "C" and 
compare it to the free fall motion, and to a motion of a fictional avalanche that started at 
some height in the building, say, H* and continued for H\. 

Here we recall that a difference between the avalanche and the free fall is that in the 
avalanche it is the top section that destroys a part of the base in its path and then adsorbs 
it, while in the free fall the base is quickly converted into the free floating chunks so the entire 



buildin 

Seffen 



g simultaneously free falls. Also, as discussed earlier, Bazant and Verdure 



, and 



9[ proposed a highly speculative mechanisms which net-result is that the avalanche 
by the top section feels almost-zero resistance when crushing, presumably, intact structure 
of the base in its path. NIST investigators imply the same in the quotation we stated earlier. 

We start by stating the equation of motion of a 0-opposition avalanche that at time t has 
dropped to position Z, 

i(zz)=gZ, (21) 



dt 

where for simplicity we neglect the effects of compaction k. Assuming that the motion at 
t = starts from rest, Z(0) = 0, at drop Z(0) = Z , Eq. fT2~Tj) can be integrated once, 
yielding a relationship between the current velocity and the position. The time t it takes 
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FIG. 6: (Phase I) Trajectories for the theoretical models, free fall (green, solid line for position, 
dots for acceleration), an 0-opposition avalanche from Hi = 28 m to the ground level (red), an 
0-opposition avalanche from H2 = 68 m (orange) and the actual descent curve from the data set 
"C" (black dots). Also given is the SAE between the descent curve and a model. In terms of SAE, 
the free fall is the best fit to the descent curve. SAE over Phase I for either avalanche is greater 
than that for the three-zone model for an entire descent (~ 3.84 m). Furthermore, discrepancies 
between the free fall and the descent curve are the greatest at the first 1.5 s, during which the 
descent data is known to be imprecise due to a low resolution of the recorded video. Conversely, 
the avalanches show systematic departure from the descent curve in the second half of Phase I, 
where the descent curve is (relatively) more precise. 




the avalanche to propagate from Zq to Z\ is given by 



= - L82 V7 + 2 ' 45 V7- 2Fl V 6'2'6'Z?" 



1 1 5 Zl 



(22) 



where = 2-P1 ( a > b; c; z) is the hyper-geometric function llj . Similarly, the acceleration 
as a function of position for Z e [Z Q , Zi] is given by 



Z= 9 - fl+2§f 



(23) 
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Finally, we use Eq. (122]) and Eq. (1231) to find the drop of the avalanche and its acceleration at 
the times at which data set "C" was recorded. For completeness we consider two avalanches, 
the first, that starts at Z\ = H — Hi and propagates to the ground level, Z2 = H, and the 
second, that starts at Z\ = H — if 2 and propagates for H\ to Z2 = H — if 2 + H\. 

In Fig. [6] we show the data set "C," the respective theoretical trajectories, and their 
accelerations and SAE. Here, SAE is calculated over the first 13 points of the data set 
"C" and the positions predicted by each theoretical model. We note that the free fall has 
SAE~ 2.09 m, on one hand side, and on the other, that SAE for the entire three-zone crush- 
up model over 25 points (~ 3.82 m) is less than SAE for either of the 0-opposition avalanche 
models (~ 4.69,4.94 m). This finding reaffirms our previous conclusion that Phase I is a 
free fall for H\ ~ 28 m and not an avalanche that started somewhere in the building and 
propagated for the same distance. 

In light of our discussion of the secondary and the primary zone we speculate what must 
happen to the base for a free fall of the top section to be possible: 

1. The core of the base is destroyed in the same sequence as the core of the primary 
zone. The NIST investigators appear to be hypothesizing this to be the case: what 
they believe is a cascade of floor failures may in fact be a staged destruction (severing) 
of the core columns in the primary zone and then continues throughout the base. 
However, the damage to the base is more extensive than the damage to the primary 
zone in that in the base the floors and their web of trusses are destroyed as well. The 
top section is later, at t — 0, released by severing the perimeter columns at Hi. These 
columns offer little to no resistance to the falling top section due to their marginal 
position and small cross section. 

2. The destruction of the entire base at t — is inconsequential to an earlier destruction 
of the core columns in the primary zone. The strength, of otherwise intact, base is 
sufficient to arrest the fallout of staged destruction of the core columns in the primary 
zone. Here, the base being annihilated is what releases the top section. 

We note that in terms of the apparent weight W the building exerts on the Earth's crust 
during its collapse, Eq. ([TBI , the two cases differ. There are two type of terms contributing 
to the apparent weight: "arrest," created by a large chunk of the building coming to a stop 
after hitting the ground or a part of the building in its path, and "release," created by a 
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large chunk breaking off the building and starting a free fall. It is a reasonable assumption 
that the seismic signal is excited by changes in the apparent weight of the building, 5w, 
given by 5w = f AW/(M g), where / is the sampling rate, while A is a difference operator 
acting on a time series of W collected at the sampling rate. 

It can be shown that if the interior, core and floors, of the base and the core of the primary 
zone are destroyed prior to the release of the top section (Case 1) than the peak in Sw from 
the first release (the top section being allowed to free fall) is comparable to, possibly weaker 
than, the peak of the first arrest (the top section reaching the ground). On the other hand, 
if the destruction of the entire base marks the release of the top section (Case 2) then the 
peak of the first release is much stronger than the peak of the first arrest. 

We believe it is the seismic signal of the collapse that can be used to deduce which of the 
two cases is more likely to have had occurred. We leave this analysis to future publications 
with our collaborators. Given our current knowledge [ljj], we favor Case 1. 
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IV. CONCLUSION 

N I 

| Secondary Zone 
| Primary Zone 
| Base 

FIG. 7: Not-to-scale illustration of the first few phases of collapse of WTC 7. In the initial building, 
left panel, during Preparatory or Null Phase the core columns are destroyed between the heights 
H 2 = 68 m and H ~ 28 m, middle panel, which splits the top section in two: the 60% damaged 
primary, and the intact secondary zone. The interior of the base is destroyed, as well. Phase I, 
right panel, starts when the perimeter columns in the base are severed at Hi, allowing the top 
section to free fall. Phase II, the "crush-up," begins when the top section reaches the ground (not 
shown) . 

Given the descent curves, the results of physical modeling, some video evidence and 
the damage review by the NIST investigators we conclude that the collapse of WTC 7 is 
comprised of four phases: 

• Phase N: Null or Preparatory phase starts 8, or so, seconds before the collapse. 
During that phase, we argue, the core between H 2 — 68 m (15th floor) and H x ~ 28 m 
is destroyed together with the base interior. 

The appearance of the building during that period, which features, among others, 
sinking of the penthouses on the top into the building, is consistent with severing of 
the core columns below H 2 . The sinking results from the sections of the core columns 
above H 2 being left suspended from the hat truss and the perimeter columns. That 
these hanging sections of core columns in the secondary zone are not destroyed becomes 
apparent during Phase III when the top section in its last moments regains its full 
(local) strength. 

• Phase I: Free Fall phase begins at t — with a sudden and total annihilation of 
the base (part of the building between the ground level and Hi). This allows the top 
section (part of the building above Hi) to free fall to the ground. 
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• Phase II: "Crush-up" begins t ~ 2.3 s into the collapse when the top section reaches 
the ground. For the next ~42 m the primary zone of the top section, which was 
compromised during Phase N, is destroyed in collision with the ground. 

• Phase III: "Crush-up" of the top section continues for the next t ~ 3.8 s as the 
secondary zone is being destroyed. While the top section now begins to decelerate, 
this, in itself, is not sufficient to arrest the collapse. The phase continues some 7.8-8.8 s 
into the collapse when the last remains of the building fall on the ground. 

We conclude that the building was destroyed in a highly controlled fashion and, contrary 
to the common sentiment, did not spontaneously collapse. 
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APPENDIX A: DATA SET "C 



TABLE II: Raw data comprising the descent curve "C." [3] It measures a height of a point on the 
building, which exact position is identified from the NIST report. 
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"We use value 155.80 m instead, and posit that the uncertainty in the drop is ±0.2 m. Only effect of this 
action is that SAE is reduced by 0.2 m by hand. 



